Trilinear Gauge Boson Couplings in the MSSM 



E. N. Argyres a *, A. B. Lahanas ^ , C. G. Papadopoulos cd ** 

and V. C. Spanos c * 

a NRCPS, 'Democritos', Aghia Paraskevi, GR 15310 Athens, Greece 
b University of Athens, Physics Department, Nuclear and Particle Physics Section, 

GR - 15771 Athens, Greece 
c University of Durham, Physics Department, Science Laboratories, 
South Road, Durham DH1 3LE, England 
d CERN, Theory Division, CH-1211, Geneva 23, Switzerland 



Abstract 

We study the C and P even WW^ and WWZ trilinear gauge boson vertices (TGV's), 
in the context of the MSSM assuming that the external W's are on their mass shell. 
We find that for energies ^Js = \J~cj 1 < 200 GeV squark and slepton contributions to 
the aforementioned couplings are two orders of magnitude smaller than those of the 
Standard Model (SM). In the same energy range the bulk of the supersymmetric Higgs 
corrections to the TGV's is due to the lightest neutral Higgs, ho, whose contribution 
is like that of a Standard Model Higgs of the same mass. The contributions of the 
Neutralinos and Charginos are sensitive to the input value for the soft gaugino mass Mi/ 2 , 
being more pronounced for values M1/2 < 100 GeV. In this case and in the unphysical 
region, < i/i < 2M\y, their contributions are substantially enhanced resulting in large 
corrections to the static quantities of the W boson. However, such an enhancement is 
not observed in the physical region. In general for 2My/ < a/s < 200 GeV the MSSM 
predictions differ from those of the SM but they are of the same order of magnitude. To 
be detectable deviations from the SM require sensitivities reaching the per mille level and 
hence unlikely to be observed at LEP200. For higher energies SM and MSSM predictions 
exhibit a fast fall off behaviour, in accord with unitarity requirements, getting smaller, in 
most cases, by almost an order of magnitude already at energies yfs ~ 0.5 TeV. 
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1. Introduction 



The Standard Model (SM) has been remarkably successful in describing particle inter- 
actions at energies around ~ 100 GeV. Precise measurementsat LEP provided accurate 
tests of the standard theory of electroweak interactions , but we are still lacking 
a direct experimental confirmation of the non-abelian structure of the standard theory. 
The WWj, WWZ and ZZ'j couplings are uniquely determined within the context of 
the SM and such couplings will be probed in the near future with high accuracy The 
study of the trilinear gauge bosons vertices (TGV's) in the e + e~ — > W + W~ process is 
the primary motivation for the upgrading of LEP200 and the potential for measuring 
these has been discussed in detail BB 0. So far there are no stringent experimental 
bounds on these couplings @ and the efforts of the various experimental groups towards 
this direction are still continuing. Precise measurements of these vertices are of vital 
importance not only for the SM itself but also for probing new physics which may be 
revealed at scales larger than the Fermi scale. 

The gauge boson vertex has been the subject of an intense theoretical study during 
the last years. In particular the WWV vertex (V = 7 or Z) has been analysed in 
detail within the framework of the standard theory as well as extensions of it, and 
its phenomenology has been discussed. The lagrangian density describing the WWV 
interaction is given by 

C wwv = -igwwvWlW^V" - h.c.) + k v W^W v F^ + ^-W^WJfV^ + ... ]. (1) 

In this gwwv = e for V = 7 and gwwv = e cot 6w for V = Z. The ellipsis stand for P 
or C odd terms and higher dimensional operators. In Eq. (1) the scalar components for 
all gauge bosons involved have been omitted, that is d-W = d-V = 0, since essentially 
they couple to massless fermions Q. At the tree level Ky and Ay have the values Ky — 1, 
Ay = 0. However radiative corrections modify these, the order of magnitude of these 
corrections being O(-) ~ 10~ 3 . Sensitivity limits of this order of magnitude will not 
be reached at LEP200 but can be achieved in future colliders where the TGV's can be 
studied in detail and yield valuable information not only for the self consistency of the 
SM but also for probing underlying new physics. Any new dynamics whose onset lies in 
the TeV range modifies Ky, Ay and deviations from the SM predictions are expected. 

Supersymmetry (SUSY) is an extension of the SM which is theoretically motivated 
but without any experimental confirmation @. The only experimental hint for its ex- 
istence derives from the fact that the gauge couplings unify at energies ~ 10 16 GeV, if 
we adopt a supersymmetric extension of the SM in which SUSY is broken at energies 
MsusY^Oil) TeV ® . Supersymmetric particles with such large masses can be pro- 
duced in the laboratory provided we have very high energies and luminocities. However, 
their existence affects Ky and Ay, even at energies lower than the SUSY production 
threshold making them deviate from the SM predictions. Therefore the study of these 

1 For on shell W's we have (□ + my^W^ = and certainly d-W = 0, while for the photon d-A 
vanishes on account of current conservation. 
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quantities may furnish us with a good laboratory to look for signals of supersymmetry at 
energies below the SUSY production threshold. In any case such studies serve as a com- 
plementary test, along with other efforts, towards searching for signals of new physics 
and supersymmetry is among the prominent candidates. 

In the SM Ky and Ay as functions of the momentum q 2 carried by the V boson 
(V — 7, Z), for on shell W's, have been studied in detail but a similar analysis has 
not been carried out within the context of the MSSM. Only the quantities K 7 (g 2 = 0), 
A 7 (g 2 = 0) have been considered @, which are actually related to the static magnetic 
dipole (fiw) and electric quadrupole (Qw) moments of the W boson. To be of relevance 
for future collider experiments the form factors k 7i z, A 7i z should be evaluated in the re- 
gion q 2 > 4M^/. The behavior of A 7i z in this physical region may be different from 
that at q 2 = 0, especially when the energy gets closer to Msusy and supersymmetric 
particles may yield sizeable effects due to the fact that we are approaching their thresh- 
olds. In those cases an enhancement of their corresponding contributions is expected, 
unlike SM contributions which are suppressed in this high energy regime. 

In this work we undertake this problem and study the C and P even WW^f, WW Z 
vertices in the context of the MSSM, with radiative electroweak symmetry breaking, 
when the external W bosons are on their mass shell. Such studies are important in view 
of forthcoming experiments at LEP200 and other future collider experiments which will 
probe the structure of the gauge boson couplings and test with high accuracy the predic- 
tions of the SM. If deviations from the SM predictions are observed these experiments 
will signal the presence of new underlying dynamics at scales larger than the Fermi scale. 



1. The SM contribution to Ak v (Q 2 ), AQ V (Q< 

Alth oug h the SM contributions to the TGV's have already been calculated in the litera- 
ture EJJ, for reasons of completeness we shall briefly discuss them in this section paying 
special attention to the contributions of fermions and gauge bosons. 

In momentum space the most general WWV vertex (V = 7 or Z) with the two Ws 
on shell and keeping only the transvers degrees of freedom for the 7 or Z can be writen 
as0 

AO Q 2 
+2 Ak v (g ail Qp - g^Q a ) + 4 A„(Q a Q p - ^-g a p)} + ••• (2) 

where Akv = K v + Ay — 1 , AQ V = —2 Ay. The kinematics of the vertex is shown in 
Fig. 1 . 

The ordinary matter fermion contributions at Q 2 = have been studied elsewhere 
0' . However in those works there is an i mpo rtant sign error which affects sub- 
stantially the results given in those references This has also been pointed out 



independently in ref. [15j . The consequences of this for the static quantities of the 
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W boson nw, Qw has been discussed in detail in ref. fll2]. The fermion contributions 



to Ak^(Q = 0) are large and negative partially cancelling contributions of the gauge 
bosons and standard model Higgs which are positive. 

At Q 2 7^ and in units of g 2 /l6n 2 the fermion contributions of the graph shown in 
Fig. 1, when the internal lines are (Pi, P 2 , P 3 ) = (/, /, /'), are given by, 

Ak v = -c v T(C g / dt / da{g f L [t 4 + t 3 (-l + R« - R f ) + t 2 (R f - R«) 
Jo Jo 

AO 2 1 
- t 3 (3t-2)a(l-a)}+g f R [R f t 2 }}-j (3) 



AQ y = -c v TiC g g J L / / da — 2 (4) 

JO Jo Lj 

where c 7 = 1, c z = R with R = M 2 Z /M^. In Eqs. <g) and (§ 

L) = t 2 (l - -r^-a(l - a)) + t(P/ - P/> - 1) + P f + ie , 

2 

and C 9 is the color factor (1 for leptons, 3 for quarks). Throughout this paper P, = jA- 

w 

where the index % refers to the relevant particle in each case. The couplings g L R appearing 
in Eqs. (|) and (|) are g{ = g s R = Q{ m for V = 7, and g[ = Q f wL , g f R = Q f wR for V = Z. 
Q{ m are the electromagnetic charges, and Q^ wLR are the weak charges for left/right 
handed fermions defined by the relation =T( — Q( m sin 2 9w- P3 is the weak isospin 
of the fermion /, ie. T( = —1/2 for the left handed charged leptons etc. 

Regarding the contributions of the gauge bosons to Akv, AQv ,(Q 2 7^ 0), one faces 
the problem of gauge dependence. This is reflected in the loss of the crucial properties 
of infrared finiteness and perturbative unitarity. Analytic calculations of Akv, AQv 
were carried out in ref. |TT| in the 't Hooft - Feynman gauge. In that reference it was 



shown that Akv h as bad high energy behaviour, growing logarithmically as the energy 
increases, violating therefore unitarity constraints. At the same time Akv exhibits an 
infrared (IP) singularity the only one surviving among several others that cancel each 
other. 

The problem of gauge dependence is simply related to the fact that in general n-point 
functions (propagators, vertices, etc) are not gauge invariant objects; only the S'-matix 
elements have this property, being physically measurable quantities. Nevertheless there 
are several ways to extract the physically relevant information from their calculation. 
The most trivial example in our case, is to take the limit Q 2 — » 0, which defines the 
static properties of W boson. Less trivial is to consider special combinations of AA; 7 
and Akz which are directly measurable in the experiment. This is for instance the case 
of e^e R — > W + W~ reaction, where the combination A/c 7 — AQ 2 /(AQ 2 — m 2 z )Akz can 
be measured. As is rather evident this combination is free of IR and gauge- dependence 
problems, as can be easily verified from the calculation of Akv [P • 
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Of course, the gauge-dependence cancels out when the full one-loop matrix element 
for a specific process, for instance e + e~ — > Vy + Vy~, is performed. This means that the 
gauge-dependence of the three-point functions is eliminated from the corresponding con- 
tributions in the four-point functions. Therefore one can project out from the four-point 
functions the relevant tensor structures multiplying Aky and AQy and after collecting 
all possible contributions to the specific process a gauge invariant and physically relevant 
result would be obtained. Nevertheless this procedure would require the knowledge of 
the full one-loop contributions for each specific process. On the other hand one can 
proceed in this direction by employing special field theory methods such as the Pinch 
Technique (PT) in order to define gauge-invariant and process-independent quantities. 



The PT has been introduced by Cornwall l!^ and has since been employed in various 



physical problems ( see ref. [16| — |19|). The main feature of the PT algorithm, is that 
it enables us to define gauge-parameter independent self energies and vertices which 
satisfy QED-like Ward identities. More recent studies have explicitly demonstrated the 
universality of the PT-algorithm for gauge boson two point functions £B. Another ap- 
proach, through which gauge-parameterQ and process independent quantities can also be 
obtained, is the Background Field Method (BFM) 0. In principle, BFM provides a 
complementary and natural way to explain why the PT rules yield vertex functions that 
have the desirable properties and satisfy simple Ward identities. 

More specifically in our case PT provides us with an algorithm to isolate from box 
graphs pieces that have a vertex like structure and allot them to the usual vertex graphs. 
The details of this proceedure and the way one extracts the relevant pinch parts Afc 7jZ 
can be traced in ref. ]T7| . Once the pinch contribution are taken into account the gauge 
boson contributions to TGV's become gauge independent approaching zero values as Q 2 
increases, as demanded by unitarity, being also free of infrared singularities. 

Furthermore, by an explicit calculation, it can be shown that the pinch parts do 
not contribute when the incoming electron carries right handed chirality, as expected. 
In fact in that case the particular combination entering the amplitude for the process 
e t e R ~~ *■ W + W~ has vanishing pinch part and hence in such a process only the genuine 



vertex parts contribute[ll 



In the sequel and in order to keep our discussion independent of the specific process 
as much as possible we shall consider the SM contributions with the pinch parts properly 
incorporated into the usual vertex contributions. As yet, there is no systematic study 
of how well the PT approximates the coefficient of the relevant tensor structures of 
the full one-loop result in specific processes. Nevertheless, we shall use it as a basis for 
comparison of the supersymmetric contributions which are well behaved and not plagued 
by gauge dependent or infrared pathologies. 



2 With respect to gauge transformations of the background field. 
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2. The MSSM contribution to Ak v (Q 2 ), AQ V (Q 2 ) 

The MSSM is described by a Lagrangian @ 

£ = £-SUSY + £so/t (5) 

where Csusy is its supersymmetric part and C S oft the SUSY breaking part. A convenient 
set of parameters to describe low energy physics is given by 

m , A , M 1/2 , tan/3(M z ) , m t (M z ), (6) 

where m , A , Mi/ 2 are the soft SUSY breaking parameters, m t (M z ) is the value of the 
"runing" top quark mass at the Z - boson pole mass M z , and tan f3(M z ) the ratio of 
the v.e.v's of the Higgses H 2 ,Hi. Complete expressions for the RGE's of all parameters 
involved can be found in the literature and will not be repeated here 0. 

At Q 2 = the MSSM contributions to AA; 7 , AQ 7 have been studied elsewhere and 
their dependences on the soft breaking parameters A, mo, M 1 / 2 , tan/3 and top quark 
mass m t have been investigated E^l. In summary the conclusions reached in that paper 
are as follows : 

i) Squarks and Sleptons have negligible effect on the dipole and quadrupole moments as 
compared to contributions of other sectors. 

ii) The bulk of the MSSM Higgs contributions to AA; 7 , AQ 7 is due to the lightest CP - 
even neutral ho, whose one loop corrected mass does not exceed ~ 140GeV. 

iii) The Neutralinos and Charginos under certain conditions may be the dominant source 
of substantial deviations from the SM predictions. This happens for values of the soft 
gaugino mass M1/2 <C Ao,mo. In some extreme cases their contributions can even 
saturate the LEP200 sensitivity limits. 

The question of whether sizeable deviations from the SM predictions can show up 
at nonzero values of Q 2 , accessible at LEP200 or other future colliders, is our principal 
motivation to study the behaviour of aforementioned quantities in the Q 2 7^ region. 

In the MSSM gauge boson and ordinary fermion contributions are like those of the 
SM and will be treated in the way desrcibed in the previous section. The additional 
contributions from q, I (squarks, sleptons), Z, C (neutralinos, charginos) as well as the 
supersymmetric Higgs contributions to Aky, AQy can be deduced from the triangle 
graph shown in Figure 1 where Pi.,2,3 are the appropriate internal lines. We will not 
consider graphs that yield vanishing contributions to both Ak v and AQ V . 

In the following we shall consider the contributions of each sector separately. 



Squarks— Sleptons (g, /) 

We first consider the contributions of the sfermion sector of the theory which can be 
read from the diagram of Figure 1 with the assignment Pi = P% — f, P3 = /' to the 
internal lines. / and /' denote sfermions. Unlike matter fermions this graph involves 
mixing matrices due to the fact that left fi and right /_£ handed sfermions mix when 
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where 



electroweak symmetry breaks down. Such mixings are substantial in the stops, due to 
the heaviness of the top quark, resulting in large mass splittings of the corresponding 
mass eigenstates £i )2 . 

In units of g 2 /l6n 2 the sfermion contributions are given by: 
Ak v = -2C a c v T( £ (Aj/K^K^ ) {k( x kI x ) 

i,j,k=l 

,i ,i t 2 (l-t)[2t-l + Rja + R f (l- a )-R f/ ] 

x dt da ^ ^ Jjs - 7 

Jo Jo Ljj, 

AQ V = 8C g c v Ti ]T (4> J ' kIkI){kI x kI) 

i,j,k=l 

/■i /■! , t 3 (l -t)a(l -a) 
x / dt / da— " (8) 

//' 

AO 2 

L% *, = t 2 + (Rf a + R f (1 - a) - l)t + R fl (1 - t) - -^-t 2 a(l - a) + ie. 

w 

In the expressions above are given by, Afj = Q^S^—T^ K( 2 Kj 2 , Ajf = Q{ m Sij 

where Q^ w and Q^ w = T/ — Q{ m sin 2 6^ are the electromagnetic and weak charges of the 
sfermion / respectively. The remaining factors c 7 , cz , and C g were defined in the 
previous section. and f[ 2 denote the mass eigenstates while K-^' diagonalize the 
corresponding mass matrices (for notation see ref. [0). In the stop sector for instance, 

~T 

the diagonalizing matrix is defined as K*.M|K* = diagonal(m 2 ^,m 2 ^). In the absence 

of SUSY breaking effects mj z, = nifji and K^' become the unit matrices. In that 
limit AQv given above cancels against the corresponding fermionic contribution of Eq. 
(|) as it should. 

Neutralinos— Charginos (Z, C) 

The neutralino and chargino sector is perhaps the most awkward sector to deal with 
due to mixings originating from the electroweak symmetry breaking effects. The two 
charginos Ci (Dirac fermions) and the four neutralinos Z a (Majorana fermions) are 
eigenstates of the mass matrices M.q and M.-^ whose explicit expressions can be found 
in ref. [[HJ. If the U, V matrices diagonalize A4q, i.e. UAi^V^ = diagonal, and O (real 
orthogonal) diagonalizes the real symmetric neutralino mass matrix Aij^, T J\4^0 = 
diagonal, then their electromagnetic and weak currents are given by 

i a,i 
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where 

= ~^f 3a u* 2 - o 2a u:, , c L m = +^o 4a v* - o 2a v* (9) 

4 = [cos 2 9 w 5 t3 -^V t2 V* 2 5 hL + U* 2 U j2 5 hR )} , h = L,R (10) 

B a/3 = ^{0 3a O w - Aa O^) , B* /3 = -B^. (11) 

Pr^l are the right /left— handed projection operators (1 ± 7s)/2. 

The contributions of this sector to A/c 7 , AQ 7 , which are given below, are calculated 
from the graph of Figure 1, with the following assignments to the internal lines -Pi,2,3, 

(P 1 ,P 2 ,P 3 ) = (C,C,Z) : 

AA; 7 = -Y f 1 dt f 1 da{F ai [t 4 + (R a - Ri - l)t 3 + (2Ri - R a )t 2 

+^t%3t - 2)a(l - a)} + G a ^{At 2 - 2t)}±- (12) 
AQ, = -8j2F ai Cdt / W 3(1 ~^ (1 ~ Q) , (13) 
where F m =| C* | 2 + | | 2 , G m = (C&C7£ + /i.e.), and 



4Cr 

L| = t 2 + {R l -R a -l)t + R a --^ 2 -i 2 a{l-a)+te. (14) 



In the equations above the indices a = 1,2,3,4 and i = 1,2 refer to neutralino and 
chargino states respectively. Note that we have not committed ourselves to a particular 
sign convention for the masses rrii, m a appearing in the sum in the equation above for the 
A/c 7 . Chiral rotations that make these masses positive also affect the rotation matrices 
and should be taken into account. 

For the couplings Akz, AQ Z we have the following contributions 

i) (P l ,P 2 ,P 3 ) = (C,C,Z) : 

Ak z = ~R Y dt da{Sf jc \t 2 {\ - t)(-t + R i a + R j (l -a) - R a ) 

+ T7H 3 ( 3t - 2)«(1 " «)] + ^^t 2 

M 2 , Myy 

-^Wa + T 5a)W° + ? ~ ( 15 ) 

A 0z = -8B £ /' * /' da ' 3 d -'Ml-«) (16) 



7 



where 

qL(R) _ /^iL*siL aL(R) . / r _i. p\\ rpL(R) _ /s-tL*s-iR /\L(R) , / t _± r>\\ 

^ija = { L/ ai L/ aj /i ji +{^^ri)) , l ija = ai^ aj^-ji +{^^ri)) 

AO 2 

Lj ja = atRi + t(l - a)Rj + R a (l -t) -t(l - t) - j^-t 2 a(l - a) + ie 
The indices i,j refer to charginos and a to neutralino mass eigenstates. 



ii) (P 1 ,P 2 ,P 3 ) = (Z,Z,C) : 

This is the same as the previous graph with {i, j, a} replaced by {p, a, i} and S*f p R \T^ p R \ 
L 2 j a replaced by the following expressions: 



S' p L f ] = -{C L ;C L ci B L p ^ + (L # R)) , T%*> = -{CZC&W + (L # R)) 

AO 2 

L' 2 ai = atR p + t(l-a)R a + Ri(l-t)-t(l-t)--^t 2 a(l-a)+ie 



In these a, p refer to neutralinos and i to chargino mass eigenstates. 



Higgses (H ,h ,A,H ± ) 

There are five physical Higgs bosons which survive electroweak symmetry breaking. Two 
of these, H and ho, are neutral and CP even, while a third A, is neutral and CP odd. 
The remaining Higgs bosons, H±, are charged. At the tree level the lightest of these, 
namely ho, is lighter than the Z gauge boson itself. However it is well known that 
radiative corrections, which are due to the heavy top, are quite large and may push 
its mass above M z . h turns out to yield the largest contributions of all Higgses to the 
TGV's, since the remaining Higgses have large masses of the order of the SUSY breaking 
scale. At the tree level the masses of all Higgs bosons involved are given by the following 
expressions : 

m " A = ~Ja% ' m H±= m A + M w (17) 
m 2 HoM = l -{{m\ + M 2 z f ± + M 2 f - AM 2 z m 2 A eos\2(3)} (18) 

The Higgs contributions can be expressed in terms of their masses and an angle 0, which 
relates the states Si = cos (3 (Real Hf) + sin (3 (Real #2) , S 2 = — sin (3 (Real Hf) + 
cos/9 (Real if°) to the mass eigenstates h ,H . The state Si is the SM Higgs boson, 
which however is not a mass eigenstate since it mixes with S 2 - When sin 2 6 = 1 such a 
mixing does not occur and ho becomes the standard model Higgs boson Si. 
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The various contributions of the Higgs bosons to Afc 7 , AQ y , are given below. The 
assignment to the internal lines (Pi, P 2 , P3) in each case is explicitly shown: 

i) (P 1 ,P 2 ,P 3 ) = (H + ,H^,A) : 

Ak 7 = D 2 (R A ,R H+ ) , AQ 7 = Q(R A ,R H+ ) (19) 

ii) (Pi, P2, P3) = (W+, W-, h ) + (H+, H-, h ) : 

AA; 7 = sin 2 9 D 1 (R ho ) + cos 2 9 D 2 (R ho , R H+ ) (20) 

AQ 7 = sin 2 #Q(i4 ,l) + cos 2 9 Q(R ho ,R H+ ) (21) 
Hi) (Pi, P 2 , P 3 ) = (W+, W-, Hp) + (H+, H-, Hp) : 

As in ii) with Rh — > Rr and sin 2 9 # cos 2 9 (22) 

where sin 2 9 = (M\ + Mf sin 2 2/5 - Ml o )/(M 2 Ho - M|J . The functions D h2 , Q appear- 
ing above are given by 

, v If 1 ,/" 1 , 2t 4 + (-2 - r)t 3 + (4 + r) t 2 
Dx(r) = - / dt / da ^ j4 ^ 23 

W 2io io t2 + r (l_ t )_ 4^t2 a(1 _ a ) + , e V ' 

W 

„ , „ N 1 r 1 , f 1 , 2t 4 + (-3 - r + P)t 3 + (1 + r - P)t 2 . . 

D 2 (r,R) = - / dt / da - - ^4 - 24 

y ! 2 Jo Jo t 2 + (-l-r + R)t + r-j^t 2 a(l-a)+ie V ; 

W 

Q(r, R) = 2 f 1 dt f 1 da t 3 (l - t)a(l - a) 

V ' Jo Jo t 2 + (-l-r + R)t + r-j&t 2 a(l-a)+ie V ' 

W 

For the Akz, AQz form factors we get the following contributions: 
i) (P 1 ,P 2 ,P 3 ) = (W + ,W_,ho + H ) : 

Ak z = i{(cos 2 9) f 1 dt f 1 da [(4 - 2R)t 4 + (R- 2)(R Ho + 2)t 3 
4 Jo Jo ' 

+ (2R Ho - RR Ho + 8 - 2i?)t 2 ]- 2 - + (sin 2 9) x (if -> h )} (26) 



AQ Z = (2 - R) {(cos 2 9) f 1 dt f 1 da 

Jo Jo 



t 3 (l-t)a(l-a) 



+ (sin 2 #) x (Ho^h)} (27) 
where L 2 H = t 2 + i2 Ho (l - i) - |^t 2 a(l - a) + ie. 

W 

ii) (Pi,P 2 ,P 3 ) = (H + ,H_,A + Ho + h ) : 

Ak z = (?-^){D 2 (R A , R H+ ) + sin 2 D 2 (R Ho , R H+ ) + cos 2 D 2 (R ho , R H+ )} (28) 
AQz = (^-^){Q(R A , Rh + ) + sin 2 Q(R Ho , R H+ ) + cos 2 Q(R ho , R H+ )} (29) 



iii) (Pi, P 2 , P 3 ) = (Hp + h , A, H+) + (A, H + hp, H+): 



R r 1 r 1 

Ak z = — {(sin 2 #) dt da 
2 Jo Jo 



i t 2 (l - t)(l + R H+ - R A a - R Ho (l -a)- 2t) 



T 2 



+ (cos 2 9) x (H -> /i )} (30) 
- 1 fi a(l -a)t 3 (l -t) 



AQ Z = 2R{(sin 2 9) [ dt f da^ — %±± — ± + (cos 2 9) x (H -> /i )} (31) 

jo jo L# 



with L^ = - f) + + i2 Ho (l - a)i + R H+ (l - t) - ^t 2 a(l -a) + ie, 
and finally, 

iv) (Pi, P 2 , P 3 ) = (Z, Hp + hp, W+) + (Hp + hp, Z, W+) : 
R r 1 r 1 

Ak z = ir{(cos 2 e) dt da [-Qat 2 + (t 3 - t 2 )(2(t - 1) + (R - R Ho )a + R Ho ) 
2 Jo Jo 

+2(R - l)at 2 ]-^- + (sin 2 9) x (H -> ho)} (32) 
AQ Z = 2P{(cos 2 9) C dt C da ^—^-^—^- } + (sin 2 9) x (if -> ^o)} (33) 
with Li = (1 - t) 2 + Pat + R Ho t(l - a)- j£-t 2 a(l - a) + ie . 

In most of the parameter space the Higgses A, iJ ± and if turn out to be rather heavy 
with masses of the order of the SUSY breaking scale; therefore all graphs in which at 
least one of these participates are small. At the same time sin 2 9 has a value very close 
to unity. Thus the dominant Higgs contribution arises solely from graphs in which a ho 
is exchanged. This is exactly what one gets in the SM with h playing the role of the 
SM Higgs boson. 



The form factors considered so far develop also imaginary (absortive) parts which 
show up as we cross the thresholds of internal particles in the loop. These can be 
calculated using the ie prescription. Note also that the pinch parts discussed before 
have imaginary parts which should be taken into account as we do. The absortive parts 
contribute to physical processes too and therefore for a complete analysis their behaviour 
as a function of the energy should be studied. For lack of space in this paper we do not 
display analytic expressions for the imaginary parts. Their behaviour as a function of 
the variable Q 2 will be discussed in our conclusion part. 
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3. Numerical Analysis — Conclusions 



Both dispersive and absortive parts of the trilinear WWj, WW Z vertices can be cast 
in the form of single integrals of known functions of t and Q 2 , which also depend on the 
physical masses of all particles involved. In fact wherever double Jq 1 dt Jq 1 da integrations 
are encountered we first perform the Jq 1 da integrations explicitly and subsequently carry 
out the Jq 1 dt integrations numerically using special routines of the FORTRAN Library 
IMSL . The advantage of using this facility is that it leads to reliable results even in 
cases where the integrands exhibit fast growth at some points or have a rapid oscilla- 
tory behaviour. The inputs in these calculations are the variable Q 2 and the arbitrary 
parameters of the MSSM. 

With the experimental inputs M z = 91.188GW, sin 2 6 W = .232, a em (M z ) = 1/127.9 
and a s (M z ) = .115 and with given values for the arbitrary parameters tan j3(M z ), 
m t (M z ), A , itlq, Mi/2 we ran our numerical routines in order to know the mass spectrum 
and the relevant mixing parameters necessary for the evaluation of the form factors given 
in the previous sections. For the running top quark mass m t (M z ) we took values in the 
whole range from 130 GeV to 190 GeV, although small values of m t are already ruled 
out in view of the CDF and DO results @. The physical top quark masses emerging 
are slightly larger by about 3% 

As for the soft SUSY breaking parameters A , mo, Mi/ 2 we scan the three dimensional 
parameter space from ~ 100 GeV to 1 TeV. This parameter space can be divided 
into three main regions i) Aq — mo — M1/2 (SUSY breaking terms comparable), ii) 
A ~ m <C Mi/2 (the gaugino mass is the dominant source of SUSY breaking ) and iii) 
Mi/2 <S A ~ m (A ,m dominate over Mi/ 2 ) . Case ii) covers the no-scale models for 
which the preferable values are A = m = 0, while case iii) covers the light gluino case. 

Regarding the values scanned for the energy variable Q 2 we explored both the timelike 
and spacelike regions for values ranging from | Q 2 | = to | Q 2 \ = 10 5 M^. For the 
timelike case, which is of relevance for future collider experiments, this corresponds to 
values of yfs ranging from GeV to about 600 My/- Both in the spacelike and timelike 
energy region we have seen that as soon as ^ exceeds ~ few TeV the contributions of 
each sector separately becomes negligible, approaching zero as the energy increases in 
accord with unitarity requirements. 

Sample results are presented in Table 1 for values of (A , m , M t / 2 ) equal to (300, 300, 
300), (0, 0, 300) and (300, 300, 80) GeV representative of the cases i),ii) and iii) discussed 
previously. The inputs for the remaining parameters are tan/3(M z ) = 2, m t (M z ) = 
170GeU. The value of a/s in these tables are respectively 190 and 500Ge^, corresponding 
to the center of mass energies of LEP200 and NLC In the same table for comparison we 
give the SM predictions for Higgs masses equal to 50, 100 and 300 GeV. With the inputs 
given above the typical SUSY breaking scale lies somewhere between 2My/ and 0.5 TeV. 
Although many sparticle thresholds exist in this region, as for instance the lightest 
of the sleptons and squarks as well as the lightest of the neutralinos and charginos, 

3 Throughout this paper whenever we refer to MSSM predictions we mean both sparticle and particle 
contributions, gauge boson and ordinary fermion contributions inclusive. 



11 



especially when M1/2 is light, these thresholds do not result in any enhancement of the 
form factors Ak lt z, AQ 7i z- Increasing the value of the dominant SUSY breaking scale 
the supersymmetric contributions to these quantities become less important approaching 
zero values. 

We now come to discuss separately the contributions of sparticles and supersymmetric 
Higgses. The sfermions yield the smaller contributions even in cases where due to large 
electroweak mixings some of the squarks, namely one of the stops, are relatively light. 
The supersymmetric Higgses yield contributions comparable to those of the SM, provided 
the latter involves a light Higgs with mass around 100 GeV. The bulk of the Higgs 
contributions is due to the lightest CP even neutral h . Thus SUSY Higgs contributions 
are like those of the SM with h Q playing the role of the Standard Model Higgs boson. The 
neutralinos and charginos in some cases, depending on the given inputs, can accomodate 
light states. Their contributions in that case are not necessarily suppressed and are the 
principal source of deviations from the SM predictions. The contributions of this sector 
are sensitive to the input value for the soft gaugino mass Mi/ 2 , being more important 
for values Mi/ 2 < 100 GeV . For such values of the soft gaugino mass and in the 
region, < y/s < 2M W their contributions are enhanced, due to the development of an 
anomalous threshold, resulting to sizeable corrections to the magnetic dipole and electric 
quadrupole moments of the W boson @ . However such an enhancement does not occur 
in the physical region and these contributions fall rapidly to zero as we depart from the 
unphysical region to values of energies above the two W production threshold. 

The total contributions to the Afc 7i z, for some particular cases, both in the MSSM 
and SM are shown in Figure 2 for values of i/i ranging from OGeV to ITeV. The region 
from to 2My/ is unphysical since the external W's have been taken on their mass 
shell. At s = the quantities A/c 7 , AQ 7 are linearly related to the magnetic moment 
and electric quadrupole moments of the jy-boson. The MSSM predictions displayed in 
Figure 2 are for the particular choice (A , m , Mi/ 2 ) = (300, 300, 80), the most interesting 
of the three cases discussed previously, since it includes light neutralino and chargino 
states. We only show the fi > case. The negative /i case leads to similar results. For 
lack of space we have only displayed the dispersive part. The absortive part turns out 
to have qualitatively a similar behaviour and we will not discuss them any further. All 
form factors tend to zero fairly soon with increasing energy reaching their asymptotic 
values at energies a/s ~ few TeV in agreement with unitarity constraints. The first peak 
observed in the unphysical region ( ^/s < 2Myy ) in the MSSM is due to the anomalous 
threshold of the Neutralino/ Chargino sector discussed previously. In the physical region 
( y/s > 2Mw) the first peak observed is associated with the ti production threshold 
while the second, around 700 GeV, is due to the threshold of the heavy neutralino states. 
We should point out that if it were not for the Neutralino/Chargino sector the MSSM 
and SM predictions would differ little. This sector is therefore the dominant source of 
deviations from the Standard Model predictions provided it accomodates light states. 

The TGV's studied in this paper have been promoted to physical observables as 
being gauge independent satisfying at the same time the perturbative unitarity require- 
ments, using the PT algorithm. The question of how to extract physical information on 
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these vertices from measurements of physical processes at e + e~ and hadron colliders has 
been the subject of many phenomenological studies and has recently triggered numer- 
ous debates among physicists in working groups at LEP200 |22| and other workshops 
[D. Zeppenfeld in [p3"| 1- and WWZ vertices are parts of the experimentally 

observed amplitudes of the e + e~ — > W + W~ process whose SM radiative corrections 
have been studied elsewhere |[24|| . The accuracy of measuring the aforementioned ver- 



tices is improved with polarized beams. In this case, it is the particular combination 
Akj — 4Q 2 /(AQ 2 — m 2 z )Akz which enters the amplitude for the reaction e£e^ — > W + W~ 
|7], p5] , and this is gauge independent and free of infrared singularities. Moreover it is 
independent of the PT algorithm, as we have already discused. Its behaviour as a func- 
tion of the energy is displayed in Figure 3 for both MSSM and SM. The parameters are 
as in Figure 2. In the MSSM both fi > and fi < cases are shown. One notices that 
MSSM and SM predictions for this quantity are very close to one another apart from 
the peak which is due to the heavy neutralino threshold. Differences however are small, 
and to be detected requires high sensitivities. 

In our analysis we have focused our attention on the TGV's considering all SUSY 
contributions to them. For a complete one loop analysis of the e + e~ — > W + W~ process 
the SUSY box contributions should be taken into account. There is no a priori reason 
why these should be small. However these graphs involve exchanges of at least one 
selectron or sneutrino whose masses are large of the order of the supersymmetry breaking 
scale. Our previous considerations on the TGV's has shown that sleptons yield smaller 
contributions as compared to other sectors and especially that of the neutralinos and 
charginos. On these grounds we expect a small effect from the supersymmetric box 
graphs in which sleptons are exchanged. Indeed it has been showncHl that the effect 
of some supersymmetric boxes is quite small at LEP200 energies. At NLC energies 
however these become comparable to the TGV contributions. Therefore for a complete 
phenomenological study their contributions should be taken into account. The effect of 
the box graphs on the TGV form factors entering into e + e _ — > W + W~ is under study 
and the results will appear in a future publication. 

Our conclusion concerning the trilinear gauge boson vertices is that for energies 
2Myy < \fs < 200 GeV the MSSM predictions differ in general from those of the 
SM but they are of the same order of magnitude. To be detectable deviations from the 
SM requires sensitivities reaching the per mille level and hence unlikely to be observed 
at LEP200. If deviations from the SM predictions are observed at these energies it will 
be a signal of new underlying dynamics, which however will not be of supersymmetric 
nature. At higher energies SM and MSSM predictions fall rapidly to zero, due to unitar- 
ity, getting smaller, in most cases, by almost an order of magnitude already at energies 
^Js m 0.5 TeV. As a result, the task of observing deviations from the SM which are 
due to supersymmetry demands higher experimental accuracies as well as a complete 
theoretical treatment which properly takes care of box contributions. 
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tan/3 = 


2m t = 170 GeV 




A ,m ,M 1/2 


300, 300, 80 


300, 300, 300 


0,0,300 


ji > /i < 


li > ji < 


/! > jj, < 




y/s =190GeV 


AA; 7 


-1.989 -1.783 


-1.793 -1.818 


-1.812 -1.841 


AQ 7 


0.903 0.297 


0.523 0.525 


0.534 0.537 


Ak z 


-2.359 -2.065 


2.209 -2.196 


-2.208 -2.189 


AQ Z 


-0.366 -1.197 


0.342 0.360 


0.331 0.354 


SM 


Ak^ = -2.005, - 


1.735,-2.118 AQ 7 


= 0.524,0.530,0.503 


predictions 


Ak z = -1.350,- 


2.437, -1.404 AQ Z 


= 0.507,0.533,0.481 




yfi =500GeV 


Ak y 


-0.262 -0.310 


-0.151 -0.207 


-0.191 -0.259 


aq 7 


0.150 0.146 


0.030 0.041 


0.056 0.069 


Ak z 


0.115 0.238 


0.198 0.191 


0.203 0.198 


AQ Z 


0.362 0.256 


-0.407 -0.325 


-0.427 -0.352 


SM 


AA; 7 = -0.250, - 


0.168,0.046 AQ 7 


= 0.054,0.057,0.064 


predictions 


Ak z = 0.147,0.208,-0.036 AQ Z 


= 0.057,0.058,0.077 



Table 1: MSSM predictions for Ak 1; z, AQ 7i z, in units of g 2 /16n 2 , for three different 
inputs of A , m , Mi/ 2 (in GeV), at LEP2 and NLC energies. Both /i > and /i < cases 
are displayed. The SM predictions for Higgs masses 50, 100 and 300 GeV respectively 
are also displayed. 



Figure Captions 

Figure 1: The kinematics of the WWV vertex. Pi,2,3 denote internal particle lines. 

Figure 2: a) MSSM predictions for the real part of Ak 1 (solid lines) and Akz-, (dashed 
lines), in units of g 2 /l6ir 2 , as functions of the energy yHs. The inputs are (A , m , Mi/ 2 ) = 
(300, 300, 80) GeV, tan (5 = 2, m t = 170 GeV. Only the MSSM case // > is displayed. 
The vertical dotted line indicates the position of the WW production threshold, b) SM 
predictions for m t = 170GeV and SM Higgs mass, ran = lOOGeV. 

Figure 3: A/c 7 — 4Q 2 /(4Q 2 — m 2 z )Ak z as a function of the energy, for the SM (solid 
line), MSSM with /i > (dashed line) and MSSM with fj, < (dashed doted line). The 
MSSM parameters are as in Figure 2. The standard model Higgs mass is taken lOOGe^. 
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